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ABSTRACT 

Marine debris is a problem for coastal areas throughout the world, including the Gulf of Mexico. To aid the NOAA 
Marine Debris Program in monitoring marine debris dispersal and regulating marine debris practices, sea surface 
height and height anomaly data provided by the Colorado Center for Astrodynamics Research at the University of 
Colorado, Boulder, were utilized to help assess trash and other discarded items that routinely wash ashore in 
southeastern Texas, at Padre Island National Seashore. These data were generated from the NASA radar altimeter 
satellites TOPEX/Poseidon, Jason 1, and Jason 2, as well as the European altimeter satellites ERS-1, ERS-2 
(European Remote Sensing Satellite), and ENVISAT (Environmental Satellite). Sea surface temperature data from 
MODIS were used to study of the dynamics of the Loop Current. Sea surface height and MODIS data analysis were 
used to show that warm water in the core of eddies, which periodically separate from the Loop Current, can be as 
high as 30 cm above the surrounding water. These eddies are known to directly transfer marine debris to the western 
continental shelf and the elevated area of water can be tracked using satellite radar altimeter data. Additionally, 
using sea surface height, geostrophic velocity, and particle path data, foretracking and backtracking simulations 
were created. These simulation runs demonstrated that marine debris on Padre Island National Seashore may arise 
from a variety of sources, such as commercial fishing/shrimping, the oil and gas industry, recreational boaters, and 
from rivers that empty into the Gulf of Mexico. 

Keywords: marine debris, eddies, barrier island beach, Padre Island National Seashore, radar altimetry, surface 
circulation, Gulf of Mexico, NASA DEVELOP 


INTRODUCTION 

Marine debris is a problem for many of the world’s coastlines, and the Gulf of Mexico is no exception. In the 
Gulf of Mexico, clockwise-rotating areas of warm water, or eddies periodically separate from the Loop Current. 
These eddies vary in size from 100-200 km in diameter and travel westward towards the Texas coast at speeds 
ranging from 2-5 km per day. The Loop Current and the eddies directly transfer debris to the western continental 
shelf. The overall prevailing wind direction is to the west, and this force additionally drives and transports the debris 
to the western Gulf of Mexico’s inner continental shelf. From the shelf, the oblique wind directions create two 
longshore transport cells. The convergence point of these two cells is inside Padre Island National Seashore in 
southeastern Texas. Initially, the convergence point was a deposition site for submersed aquatic vegetation and 
shells. Presently, it accumulates vast quantities of anthropogenic marine debris (National Park Service, 2008). 

The beaches of Padre Island are now routinely littered with trash and other discarded items from the petroleum 
industry, commercial fishing vessels, recreational boaters, and other sources (National Park Service, 2008). This 
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debris can potentially generate a variety of human health risks to individuals who may frequent the island for 
recreational activity, such as locals or tourists. For income, Padre Island’s economy is heavily dependent on sport 
fishing, kayaking, and beach recreation, and ecotourism plays a major role in the Coastal Bend’s tourism-based 
economy (Lee, 2009). The total impact of ecotourism on household earnings is $233.5 million and on employment 
is 8,748 jobs (Lee, 2009). Marine debris on the coastline can influence the return rate of ecotourists and have a 
disastrous effect on the local economy. Padre Island, which is nearly 70 miles long, is the largest undeveloped beach 
in the world and can accumulate up to one ton of marine debris per linear mile (National Park Service, 2005b). Due 
to the island’s large extent, there is a substantial amount of marine debris accumulation, and this in turn can become 
a significant human health concern. For example, waste materials that accumulate along the seashore, like used 
medical syringes and broken glass, pose as physical hazards to humans. Other waste items, such as used condoms 
and tampon applicators, can transmit disease. Additionally, these items can affect water quality and could 
correspondingly, result in illness (Sheavly, 2007). Water that has been polluted with harmful bacteria and viruses, 
upon consumption or skin contact, “can result in infectious hepatitis, diarrhea, bacillary dysentery, skin rashes, and 
even typhoid and cholera” (Sheavly, 2007). 

The biota of Padre Island can also be threatened by the accumulation of marine debris. Padre Island is situated 
along the Central Fly way, which is an important migration route for American birds. Three hundred and eighty 
different bird species visit the island annually; this includes nearly half of all bird species documented in America, 
including migratory, overwintering, and resident bird species (National Park Service, 2011a). According to the 
National Park Service (2011a), thirteen of these species are known to be threatened or endangered. Endangerment of 
these migrating shore birds could be increased if they were to ingest marine debris that was either toxic debris or 
non-food items. Debris, upon ingestion, could abate the absorption of vital nutrients, and/or upon entanglement, can 
suffocate the birds (Ryan, 1990). Furthermore, if adult birds consumed debris and then regurgitated it to feed their 
chicks, this can lead to the death of the young and lowers nesting colony production (Laist, 1987; Fry et al., 1987). 

Throughout the Gulf of Mexico region, marine debris monitoring and cleaning up programs have been 
established, however, currents and eddies transfer marine debris over great distances, and these effects can challenge 
efficient monitoring and cleanup efforts. In 2006, former President George W. Bush signed into law the Marine 
Debris Research, Prevention, and Reduction Act (NOAA Marine Debris Program, 2011). After Hurricanes Katrina, 
Rita, and Ike, NOAA created the Marine Debris Program under former President Bush’s Marine Debris Research, 
Prevention and Reduction Act. This act established NOAA’s Marine Debris Program, among other tasks, to attempt 
to reduce the amount of marine debris in the Gulf of Mexico. Following the formation of a partnership with the US 
Coast Guard, the Marine Debris Program began to use sonar to survey the Gulf of Mexico coastline for underwater 
debris. However, the program does not currently utilize spaceborne sensors to analyze floating marine debris and its 
relationship to the Gulf of Mexico Loop Current and surrounding eddies. For large-scale monitoring of surface 
circulation processes which transport marine debris the Marine Debris program and others like it would benefit from 
the capabilities that NASA satellites such as TOPEX/Poseidon, Jason- 1, Jason-2 and other remote sensing platforms 
with radar altimeters provide. 

For this project we partnered with NOAA’s Marine Debris Program, the Padre Island National Seashore 
Education and Outreach Department, and the Colorado Center for Astrodynamics Research located at the University 
of Colorado, Boulder. The research from this project will provide NOAA’s Marine Debris Program with a better 
understanding of how the Gulf of Mexico Loop Current and surface circulation patterns affect coastal area marine 
debris dispersal. Padre Island National Seashore will use the maps, visualizations and satellite data products 
produced by this project in their education and outreach programs, such as Talking Trash. Each year the Talking 
Trash program educates over 10,000 Texas middle and high school students about marine debris and pollution by 
discussing its effects on the environment. Dr. Robert Leben and his colleagues at the Colorado Center for 
Astrodynamics Research have played a key role in developing the methodologies for marine debris monitoring by 
providing sea surface height anomaly data. These methodologies will be transferred to the project partners, enabling 
them to track and monitor the Loop Current in the future. 


STUDY AREA 

The study area for this project includes the entire Gulf of Mexico with a specific focus on Padre Island National 
Seashore, Texas (Figure 1). Padre Island is located along the southern Texas Gulf Coast (Figure 1) and “is the 
world’s longest undeveloped stretch of barrier island” (National Park Service, 2011b). The island is oriented north- 


Pecora 18 -Forty Years of Earth Observation... Understanding a Changing World 
November 14 - 17, 2011 □ Herndon, Virginia 



south, the Gulf of Mexico lies to the east and Laguna Madre lies adjacent to the west. In 1964, the artificially- 
dredged Port Mansfield Channel divided the island and as a result, the terms ’’North Padre Island” and ’’South Padre 
Island” are often used to refer to the separate portions of the island (National Park Service, 2005a). Padre Island 
includes Padre Island National Seashore. The national seashore is approximately 70 miles (110 km) long with 
65.5 miles (105.4 km) of Gulf beach. A variety of pristine beach, dune, and tidal flat environments can be found 
within the park. 



Figure 1. Gulf of Mexico region with an inset map of Padre Island, Texas. 

Study Period 

The project partners at NOAA’s Marine Debris Program expressed interest in altimeter and remotely sensed 
data from 2008-2010. Therefore, the study period for this project ranges from the beginning of 2008 to the end of 
2010. The Colorado Center for Astrodynamics Research provided SSH data for this time frame. This study period 
will also enable the analysis of the effects hurricanes, such as Gustav and Ike, have had on the dispersal of marine 
debris. After the 2010 BP oil spill incident, surface circulation processes in the Gulf of Mexico gained national 
attention. Since the oil spill occurred during the designated study period, data collected may prove additionally 
useful for studying the dispersal of oil. 


METHODOLOGY 


Data Acquisition 

Sea Surface Height Anomaly. Sea surface height anomaly data were provided by the Colorado Center for 
Astrodynamics Research (CCAR) at the University of Colorado, Boulder. This daily dataset requires the use of 
many altimeter satellites including: TOPEX/Poseidon, Jason 1, Jason 2, ERS-1, ERS-2, and Envisat. The sea 
surface height anomaly was calculated by combining all of the raw altimeter data for each satellite into separate 
cycle files. The data were referenced to a mean sea surface height model to generate sea surface height 
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anomaly. The along- track data were next interpolated to a reference ground track to insure there are always 
samples at the same geographic location for each cycle. Then, the data were detrended using a linear least 
squares smoothing filter to remove all of the low frequency data. To produce a daily product from satellites 
that have repeat cycles ranging from 10 days (Jason) to 35 days (Envisat), a specific amount of data are used 
from each satellite for the day of interest. For example, if data are being processed for January 15 th , and Jason 1 
and Envisat data are available, data is taken from Jason 1 for January 10 th -20 th and Envisat for December 28 th - 
February 1 st . A full 10-day cycle of Jason data and a full 35 day cycle of Envisat data around the day of interest 
are used. The two datasets are then combined using a weighting scheme that places less emphasis on data that 
is farther away from the day of interest. Using a multi-grid technique, the along track data is used to fully 
populate a 0.25-degree evenly spaced grid. Once the data is gridded, it is referenced to a model mean to produce 
an estimated absolute sea surface height of the Gulf of Mexico. The Colorado Center for Astrodynamics 
Research at the University of Colorado, Boulder conducted the processing to calculate sea surface height 
anomaly (Leben, 2002). 

MODIS Sea Surface Temperature. Four-kilometer spatial resolution, weekly, MODIS Sea Surface 
Temperature (SST) data was downloaded from NASA’s Oceancolor Web (NASA, 2011) as a Level 3 mapped 
product. MODIS SST data is available in a variety of spatial and temporal resolutions; spatially at 4km and 
9km, and temporally at daily, weekly (8 day), and monthly. The Level 3 mapped products used in this study are 
global gridded data sets with all points filled, even over land. The Level 3 binned products only contain bins 
with valid SST data. No land bins are kept. All mapped products are derived from the binned products. MODIS 
SST measurements are derived from the MODIS thermal infrared channels (ll-12um). MODIS SST data was 
downloaded for the time period between November 17 th , 2009 to June 17 th , 2010. 

Wind Data. Wind data used in this study was acquired from the National Data Buoy Center (NDBC)’s 
historical data archives (NOAA, 2010). The NDBC is a branch of the National Oceanic and Atmospheric 
Administration (NOAA) that, “provides high quality meteorological/environmental data in real time from 
automated observing systems that includes buoys and a Coastal-Marine Automated Network (C-MAN) in the 
open ocean and coastal zone surrounding the United States” (NOAA, 2009). Meteorological measurements 
from the buoys are recorded every day, every hour, every ten minutes (Historical wind data available at 
http://www.ndbc.noaa.gov/maps/west_gulf_hist.shtml ). 

The primary wind parameters used for this project were wind direction and wind speed. Data was acquired 
from each of three (NDBC) buoy stations-42020, 42019, & 42002-over the past five years. For buoy 42020 
the years were 2006, 2007, 2008, 2009, & 2010; for buoy 42019 the years were 2006, 2007, 2008, 2009, & 
2010; for buoy 42002 the years were 2005, 2006, 2007, 2008, & 2009 (this buoy malfunctioned during the year 
2010 ). 

Data Processing 

Geostrophic Flow Velocities. A Matlab program was written to calculate geostrophic flow velocities (Figure 2) 
using the sea surface height data that was provided by the Colorado Center for Astrodynamics Research. 
Geostrophic flow velocities were calculated as a balance of the horizontal pressure gradient and the Coriolis 
force. However, flow is also affected by other variables, such as friction. The equations used to calculate 
geostrophic flow assume that the flow has no acceleration, the horizontal velocities are much larger than the 
vertical, the only external force is gravity, and the friction is small. The equation used to calculate geostrophic 
flow follows Marshall et al. (2003): 


! g tyr i g tyr i 

Geostrophic Flow = Vu 2 + v 2 Where: ll = and V = 

f fy f & 

Where: 
g = gravity 

f — 2 Q sin(<£>) w h ere; q j s the rotation rate of the earth (7.2921 159 x 10" 5 rad/sec), and cp is the latitude 
(This equation is the Coriolis parameter) 
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h = sea surface height 
y = latitude 
x = longitude 

After this equation is run in Matlab to calculate surface velocity, the “quiver” command is used to create a 
vector field showing flow direction and magnitude. 

Particle Paths. Particle paths can be tracked by using a numerical integration method to integrate through time 
and space. The previous term DEVELOP team used Euler’s integration method to track a particle’s path by 
using the following equation (Weisstein): 


•^n+1 — ^ 

tn+1 ~ tn At 


Where x = position, v = velocity, and t = time 

For this term’s DEVELOP team, a higher-order integration method was used to more accurately track particle 
paths. The second-order Runge-Kutta method, which uses a trial step at the midpoint of an interval to cancel 
out lower-order errors, was implemented. This particular method was chosen because it provided the necessary 
accuracy needed without requiring a significant amount of processing time (Weisstein): 


Given a vector cp n of unknowns at time t n and the first order differential equation ^ = f(cp,t ), the second 
order Runge-Kutta estimate for (p n+1 is given by: 

ki ~ At * f{(p n , t n ) 
k 2 = f(<p n +-k 1 ,t n +- At) 

(pn+l = (pn + k 2 

Where At = t n+1 — t n 

In our case, cp would be calculated for both the longitude and latitude components. This equation was applied in 
Matlab and built upon the equations for geostrophic flow. 

Foretracking and Backtracking Simulations. A cluster of particles was initially placed near known oil and 
gas infrastructure. Foretracking analysis was performed to determine the possible dispersal and locations of 
these particles. The initial particles in Figure 4, (left map) are traced for a 60-day period. This procedure was 
done for each day in January 2008 and the cumulative results are shown in Figure 4 (right map). Particle 
backtracking was performed to gain a better understanding of where debris on Padre Island National Seashore 
may have originated. Below, the initial particle positions shown in Figure 5 (left map), are traced back for a 
period of 60 days. This backtracking is done for each day in January 2010 and the collective resulting particle 
positions are shown in Figure 5 (right map). 

MODIS Time Series. MODIS Sea Surface Temperature (SST) data was downloaded from NASA’s 
Oceancolor Web (NASA, 2011) as compressed archive files (.tar.gz format) and then extracted to the original 
Hierarchical Data Format (.HDF). Data was then converted from the original .HDF format into the native 
ERDAS .img format using ERDAS Imagine. Then, the data was imported into ArcMap to apply a color ramp. 
A stretched color ramp was used to symbolize the SST value. The goal of this procedure was to create an 
absolute color scale that allows direct comparison for each image. The files were then saved in .JPG format for 
later use in a time-series animation and graphics created using Windows Live Movie Maker. The animation 
consists of a simple time series slideshow of SST data images with corresponding date ranges and color scales. 

Wind Data. The meteorological program WRPLOT that was used to analyze the wind data only accepts data 
for every hour, so the raw data from the National Data Buoy Center has to be converted into an acceptable 
format (Lakes Environmental Software, 1995-2011). The WRPLOT program only accepts data within certain 
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formats, i.e. the Lakes Format, so the data was converted to this format. Since there was a large volume of raw 
data, a simple C++ program was written to automate the process of arranging the raw data files into the ‘Lakes 
Format’. For this format, the raw data is arranged by columns in the order of year, month, day, hour, minute, 
wind direction, wind speed, GDR, GST and GTIME (descriptions and measurements of this data are available 
at www.ndbc.noaa.gov/measdes.shtml). The Lakes Format required the columns be arranged in the order of 
station ID, year, month, day, hour, wind direction and wind speed (all units of measurements being the same). 
Additionally, the first line of the Lakes Format data file must state “Lakes Format”. 

After the data were converted to Lakes Format, it was imported to the WRPLOT program (Lakes 
Environmental Software, 1995-2011) which is then used to create a wind rose for visualizing the wind data. 
Then, the wind rose was exported directly to Google Earth; however, the point of location in degrees, minutes, 
and seconds must be known prior to exporting the wind rose. Wind roses for 5 -year data from each of the four 
seasons were also created. 


RESULTS 


Analysis of Results 

Radar Altimeter Data. Radar altimeter data from several NASA and European Space Agency (ESA) remote 
sensing platforms was processed to produce maps of sea surface height, sea surface height anomaly, geostrophic 
velocities, and to predict particle paths. A map of processed sea surface height and height anomaly data (Figure 
1) demonstrated that the warm waters of the Loop Current can raise the local sea level by over 30 cm in the 
location where the Loop Current flows into and exits the Gulf of Mexico. The eddies that are separate from the 
loop current are warm core eddies, and the center of these eddies also have an elevating effect on local water 
levels. The water in the core of one of these eddies can also be as high as 70 cm above the waters surrounding 
it (Figure 1). As an eddy moves across the Gulf of Mexico, it can be tracked by following the area of elevated 
water. 


SSH for 

January 1 , 2010 


30°N 
28°M 
26°M 
24°N 
22°N 
20°N 
18 °N 

98°W 96°W 94°W 92°W 90°W 88°W 86°W 84°W 82°W 80°W 


Figure 1. Example of SSH Product created from blended satellite altimeter data processed in Matlab. 

A map of geostrophic velocities (Figure 2) shows the Loop Current is a relatively fast-moving current. The 
highest velocities were mainly observed where the Loop Current enters the Gulf of Mexico through the Yucatan 
Channel (Figure 2). Relatively high velocities are also observed where the current exits the Gulf through the 
Florida Straits (Figure 2). This could be due the flow of the current becoming narrower and more channelized 
as it goes though the Yucatan Channel and Florida Straits. As the flow of the current moves farther north into 
the Gulf of Mexico, its flow field begins to expand and loses stability. This loss of stability is what ultimately 
drives the eddy shedding process. 

Pecora 18 -Forty Years of Earth Observation... Understanding a Changing World 
November 14 - 17, 2011 □ Herndon, Virginia 



Predicted Particle Path with Geostrop hie Velocity (Runge-Kutta Method) 
January 1 , 2010 - February 1 , 2010 



PU IN 

98°W 96°W 94^ 92^ 90°W SSW 86*W 84°W 32^ 80°W 


m/s 


Figure 2. Geostrophic velocity with predicted particle path using Runge-Kutta integration. 

Particle path prediction products (example shown between 26-28°N, 94-96°W, Figure 2) were initially 
produced using the Euler numeric integration technique. The Euler method proved to be fairly inaccurate due to 
the fact that some predicted particle paths left their respective eddies and re-entered the flow field of the Loop 
Current. This phenomenon is not observed in reality, and should not be occurring in the model. During the 
summer term, the DEVELOP team used a different integration technique known as the Runge-Kutta technique. 
According to advisers at the Colorado Center for Astrodynamics Research, using the Runge-Kutta method will 
improve the accuracy of the particle path prediction products. 

MODIS Sea Surface Temperature. MODIS Sea Surface Temperature (SST) data was used to analyze the 
dynamics of the Loop Current’s position, and to track its flow field and eddy shedding process. MODIS SST 
data shows the position of the loop current in the Gulf of Mexico is widely variable. The Loop Current will 
sometimes travel in a nearly direct path through the Yucatan Channel to the Florida Straits. This was seen in the 
MODIS SST data that was collected during late November and into mid December 2009. By mid to late 
December 2009 through March 20 1 0, the flow of the Loop Current had taken a more northerly path and showed 
an extension of its flow field farther north and into the central Gulf. By mid- April 2010, MODIS SST data 
(Figure 3) showed the flow of the Loop Current began to lose its stability and shed a large eddy. By the middle 
of June 2010, the eddy appeared to completely separate from the main flow path. The eddy could not be 
tracked during the summer months of 2010 using MODIS SST data because the warm waters of Gulf of Mexico 
obscured the warm waters of the Loop Current and its associated eddies. 



Figure 3. MODIS Sea Surface Temperature in Gulf of Mexico. 
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Figure 4. Image on left shows initial placement of particles before 60-day foretracking is 
performed. Image on right shows cumulative particle locations after the simulation is ran using 
each day in October, 2010, as the starting date. 



Figure 5. Image on left shows initial placement of particles before 60-day backtracking is 
performed. Image on right shows cumulative particle locations after the simulation is ran using 
each day in January, 2010, as the starting date. 

Wind Data. The time period used to create the wind roses assumed all data for the past 5 years. Wind roses 
(Figure 6) were created for three separate buoys (42020, 42019, & 42002), using annual and seasonal data. The 
wind roses each use 36 wind directions and a contrasting color scheme. The rose can then be exported to a .kml 
file and imported into Google Earth. The maps show all wind directions and wind speeds for each location. 
Maps were also created to depict a dominant seasonal wind. Although winds are variable, a dominant South 
East (SE) to South-South East (SSE) trend was observed in the five-year mean wind data (Figure 6). 

Errors and Uncertainty 

Particle Path Prediction Products. During the DEVELOP team spring term of this project, Euler’s integration 
method was used. This method has a truncation error of 0(/i 2 ). After noticing unrealistic behavior in some of 
the predicted paths, our advisors at the Colorado Center for Astrodynamics Research suggested using a higher- 
order integration method to improve the accuracy. During the DEVELOP team summer term, a higher-order 
integration method was implemented, namely the second-order Runge-Kutta method. The truncation error 
associated with this method is 0(/i 3 ). This method used the midpoint of an interval as a trial step, and as a 
result generated much smoother paths and more accurate results. 
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Figure 6. Five-year mean wind data from NDBC, processed using WR Plot software, exported as a .kml file 
and displayed using Google Earth. 

Wind Data. The C++ program accepts the filename of the raw data file as well as the station ID for the 
observed buoy as user input. The program then proceeds to arrange the data into the Lakes Format as 
previously described. After arranging the data, it then writes the output to a new file. However, there are a 
several nuances that are noteworthy. The raw data from NDBC breaks down the hour further into ten-minute 
intervals. As a result, those values were averaged to obtain a desired single value per hour. Unfortunately, the 
Lakes Format only accepts wind direction and wind speed as integer values. The program for this project 
therefore rounds the resultant averaged numbers to the nearest integer. While the overall effect this has is likely 
minimal, it should be noted that there is some amount of error involved in the final Lakes Format output. 


CONCLUSIONS 

Marine debris may negatively impact both human and marine habitats by affecting the aesthetic value of the 
beaches, creating water quality issues, and leading to the death of marine organisms. Debris on Padre Island 
National Seashore may come from many different sources, such as commercial fishing/shrimping, the oil and gas 
industry, recreational boaters, and from rivers that empty into the Gulf of Mexico. Regardless of the source, marine 
debris is an environmental issue and hazard that can directly affect coastal ecosystems and inhabitants. 
Understanding surface circulation in the Gulf of Mexico will help provide decision makers and communities 
additional information that can be used for the in elimination and/or control of marine debris. 

During the spring term, the DEVELOP team designed a methodology that was applied for the summer 2011 term. 
This methodology included the utilization of satellite altimeter and MODIS data, as well as various in-situ data. The 
methodology will help NOAA’s Marine Debris Program and the Padre Island National Park Service improve upon 
marine debris monitoring, cleanup, and education efforts. 

RADAR altimeter data provided a useful dataset for creating visual geospatial depictions of sea surface height and 
height anomaly, calculating geostrophic velocity, and predicting particle paths. The MODIS Sea Surface 
Temperature (SST) data proved useful for monitoring the Gulf Loop Current; however, during the summer months 
of 2010 the warm waters of Gulf of Mexico obscured the Loop Current and its associated eddies. Because sea 
surface temperatures in the Gulf of Mexico are universally warmer in summer months, MODIS SST data is only 
practical for monitoring the Loop Current in cooler months. Analysis of wind direction and wind speed data prove 
beneficial for understanding how marine debris is ultimately dispersed to the shoreline after eddies dissipate as they 
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interact with the continental shelf. Although wind is not the only force acting upon particles in a near-shore setting, 
it does play a key role in marine debris dispersal. 
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Improve spatial-temporal understanding of the Loop Current and its 
associated eddies, as well as their effects on marine debris 
trajectories and dispersal to coastal areas 


• Analyze debris accumulation patterns on Padre Island, Texas 

• Understand which critical marine habitats may be affected by marine 
debris pollution 

• Demonstrate the usefulness of satellite altimeter data and MODIS 
SST data in monitoring mesoscale marine surface circulation 
processes 
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TOPEX/Poseidon Jason 1 & 2 
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METHODOLOGY 


CONCLUSIONS 


Satellite altimeter data from multiple sources enabled geospatial 
maps of sea surface height, geostrophic velocity, and particle path 
predictions to be produced. 

MODIS SST data provided a means to monitor the position of the 
loop current; however, this method is not useful in the warmer 
summer months. 

Particle path predictions on surface waters can be accurately 
represented using the Runge-Kutta integration method and tri-cubic 
interpolation. 

Wind data from NDBC buoys provided useful information for 
inferring the dispersal of marine debris as it nears the coast. 

Marine Debris may negatively impact critical habitat for many 
endangered and threatened species, such as the Kemp’s ridley sea 
turtle. 


ABSTRACT 


The goal of this project was to apply satellite data to model surface circulation in the Gulf of Mexico as a means to aid in the monitoring of marine debris 
dispersal and the regulation of marine debris practices. Marine debris is a persistent problem for coastal areas throughout the world. In the Gulf of 
Mexico, the Loop Current flows north through the Yucatan Channel, loops east, then south, exiting through the Florida Straits. Clockwise-rotating areas 
of warm water, known as eddies, periodically separate from the Loop Current. These eddies have the potential to trap and transport debris onto shores, 
such as the Padre Island National Seashore. The latter is a 68 mile long barrier island beach in southeastern Texas and is the longest undeveloped 
beach in the world. This naturally formed beach can accumulate up to one ton of marine debris per linear mile. This project used sea surface height and 
height anomaly data based on NASA RADAR altimeter satellites TOPEX/Poseidon, Jason 1, and Jason 2; as well as European RADAR altimeters 
onboard ERS-1 , ERS-2 (European Remote Sensing Satellite), and ENVISAT (Environmental Satellite). This project also employed MODIS sea surface 
temperature (SST) data to aid in monitoring and visualizing the Loop Current. The sea surface height and height anomaly data were processed to 
calculate geostrophic flow velocities and predict particle paths. This research provided NOAA’s Marine Debris Program and the Padre Island National 
Seashore with a better understanding of how the Loop Current and surface circulation patterns disperse marine debris to the region. A methodology to 
monitor Gulf of Mexico surface circulation and predict particle paths by using RADAR altimeter data were provided to partnering agencies. The project 
provided maps of debris trajectories and geostrophic currents that demonstrate the usefulness of satellite altimetry data to monitor oceanographic 
processes which in turn affect the distribution of marine debris. 


RESULTS CONTINUED 


A 1994 - 1995 study conducted by 
the Minerals Management Service 
and National Park Service collected 
and catalogued 63,410 debris items 
from PINS. The commercial 
fishing/shrimping industry and the 
oil and gas industry were identified 
as the two largest contributors to 
marine debris on PINS. 


RESULTS 



Produced from June 10, 2011 
Landsat imagery 


http://develop.larc.nasa.gov 




months, warmer Gulf temperatures prevent the differentiation between 
the warm Loop Current and ambient Gulf waters. 


Because wind direction is a contributing factor to near-shore marine debris 
dispersal, wind roses are useful in understanding wind patterns and their 
potential effect on marine debris trajectories. Data source: NOAA NDBC 



Pradlclnd fttriicln Palh wilh Gootfrophic Vnlncily 
NLi* l.rflOH. Juih' 15. MOB 


Image output from Matlab showing predicted particle path 
(yellow line) overlain on geostrophic velocities base map. 
Produced from blended RADAR altimeter data. 





Same particle path prediction exported from Matlab as a 
shapefile and displayed over oil rig location and bathymetry 
map using ESRI ArcGIS 
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Community Concerns 


• Ecotourism is the fastest growing segment of the 
Texas Gulf Coast’s tourism-based economy. 

Marine debris can negatively affect the ecotourism 
industry, in turn affecting the Texas economy (Lee, 
2009). 

• Within the Coastal Bend Area of Texas, ecotourism 
employs 8,748 citizens and provides $233.5 M in 
household earnings (Lee, 2009). 

• Padre Island National Seashore (PINS) is a critical 
part of the Central Flyway, a major migratory route 
for 380 species of migratory birds (13 threatened 
or endangered) (www.nps.gov). Marine debris 
threatens this vital avian habitat. 
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Community Concerns 


• Marine debris also threatens habitat of the highly endangered Kemp’s 
ridley sea turtle. PINS is one of only two remaining global nesting sites 
for this species. Between 1996 and 2010, PINS’ Kemp’s ridley nests grew 
from 6 to 140 confirmed nests, (www.nps.gov) 


Number of Kemp’s Ridley Sea Turtle Nests Documented on the Texas Coast Since 1996 
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Partners / Collaborators 
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William Botts 

Education Coordinator - PINS 
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Partner Needs 


* Visual geospatial animations to Loop Current 
convey sources and Altimetry 

mechanisms of dispersal to 
enhance marine debris 
regulation, enforcement and 
education 




An improved method to 
monitor the Gulf of Mexic 
Loop Current and marine 
debris trajectories as well as 
predict debris particle transport 
paths 
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Study Area 



* Study Area - Gulf of Mexico & Padre Island, Texas 
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Study Area Ecosystems 



Corpus Christi 


Padre Island, Texas 
Study Area 
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Left: Padre Island General Land Cover Map 
Produced from June 10, 2010 Landsat Imagery 


% 



Above: Padre Island National Seashore Ecosystems 
Photographs Credit: National Park Service 
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Potential Sources of Marine Debris 
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Rig Locations and Major Shipping Lanes 
In the Northern Gulf of Mexico 
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A 1994-1995 study 
conducted by the 
Minerals Management 
Service and National 
Park Service collected 
and cataloged 63,410 
debris items from 
PINS. The commercial 
fishing/shrimping 
industry and the oil 
and gas industry were 
identified as the two 
largest contributors. 
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Project Methodology 





• Sea surface temperature data utilized: 

- MODIS onboard Aqua and Terra 

- Subset and converted from HDF format 
in ERDAS Imagine 

- Color ramp added in ArcGIS 


* Satellite altimeter data utilized: 

Acquired from Colorado Center for Astrodynamics 
Research, Used to Calculate Daily Sea Surface Height 

- TOPEX/Poseidon 

- Jason 1 

- Jason 2 

- ERS 1 

- Envisat 




Sea Surface Temperature (degrees Celcius) 


#MATLAB 
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Project Methodology 


Data acquired from the Colorado Center 
for Astrodynamics Research 
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Sea Surface Height (SSH) Calculated Using 
Blended Satellite Altimeter Data 
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Altimeter Data 
Processed 
in Matlab 


Geostrophic Velocity for: 
O1-Jan-201O 
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Calculated as a Balance of the Horizontal 
Pressure Gradient and the Coriolis Effect 


Particle Path Calculated by 
Integrating Over Time and Space 
Using the Runge-Kutta 
Integration Method 


Predicted Particle Path with Geostrophic Velocity (Runge-Kutta Method) 
January 1, 2010 - February 1, 2010 
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Project Methodology 


Updated Matlab code from the Colorado Center for Astrodynamics Research outputs a series of images 
showing the particle’s movement per day with sea surface height, serving as a basis for our tracking system. 




Modifications were made to output the particle 
paths into a single image with geostrophic velocity 
as well as a shapefile to be displayed using ArcGIS 
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Predicted Particle Path with Geostrophic Velocity 
May 1,2008- June 15, 2008 


Predicted Particle Path 1 r 

May 1 , 2008 - June 15 T 2008 
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Project Methodology 


Comparison of the Implementation of 
Euler and Runge-Kutta Integration Methods 



Particle Path Prediction Using the Euler Integration 
Method for January 15, 2010. In this case, the 
prediction is not an accurate representation of 
what typically happens in reality. 


Particle Path Calculated by Using the 
Runge-Kutta Integration Method for 
January 15, 2010. This is a more accurate 
prediction. 
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Project Methodology 


Data acquired from the National Data 
Buoy Center-Historical Marine Data 
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WRPLOT Program 
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Data processed 
using C++ 

into readable 
Lakes Format 



Data exported 


for analysis 


A wind rose was created using the WRPLOT program. 



Wind data analyzed for near-shore locations 


Lakes Format 


Displayed Using Google Earth 
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Results - MODIS SST 


• Because the Loop 
Current is a warm 
water current, it can 
be monitored using 
SST data. 


•During the summer 
months, this method 
has limited utility in 
monitoring the Loop 
Current. 



25 November -2 December, 2010 
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Results - MODIS SST 
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Results- RADAR Altimetry 






Image output from Matlab 
showing predicted particle 
path overlain on geostrophic 


velocites base map 


Predicted Particle Path with Geostrophic Velocity 
May 1,2009- July 1,2009 


iS 

Same particle path prediction 
exported from Matlab as a shapefile 
and displayed over oil rig location and 
bathymetry map using ESRI ArcGIS 


Predicted Particle Path 
May 1,2009 -July 1,2009 
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Results- RADAR Altimetry 
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Figure 1: 404 particles were initially placed 
near the oil and gas fields. 


Figure 2: Resulting particle positions after 
60 days foretracking 


A cluster of particles were initially placed near known oil and gas infrastructure. 
Foretracking analysis was performed to determine the possible dispersal and 
locations of these particles. The initial particles in Figure 1 were traced for a 
60 day period. This procedure was done for each day in January, 2008, and 
the cumulative results are shown in Figure 2. 
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Results- RADAR Altimetry 
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Particle backtracking was performed in order to gain a better understanding of 
where debris on PINS may have originated. Below, the initial particle positions 
shown in Figure 1 were traced back for a period of 60 days. This backtracking 
was done for each day in the month of January 2010, and the collective 
resulting particle positions are shown in Figure 2. 


Figure 1 : 1276 particles were initially placed 
near Padre Island, Texas. 


Figure 2: Resulting particle positions from 
backtracking analysis 
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Backtracking Analysis 
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Results - Wind Data 


• This image depicts annual 
winds for each buoy, over the 
last 5 years of available data. 

• Winds are shown in the 
“blowing from” direction. For 
example, buoys 42020 & 42019 
show a dominant SSE wind, 
while buoy 42002 shows a 
dominant SE wind. 

• Near-shore wind direction is 
a contributing factor affecting 
how debris is transported 
once an eddy reaches the 
continental shelf and has 
dissipated. 
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Results - Wind Data 
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During spring months, buoys 42020 & 
42019 have a dominant SSE wind, while 
buoy 42002 shows a dominant SE wind. 
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Winds during the summer months are 
slightly more variable, although S and SE 
winds remain dominant. 
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Results - Wind 



Notice in the fall months, a strong NE & 
N wind presence, as well as SE and 
SSE winds. 
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Although the dominant wind stays SE, 
a northerly wind is present during the 
winter months. 
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Discussion and Conclusion 


Conclusions: 

• Use of satellite altimeter data from multiple sources enabled the production of 
geospatial maps of SSH, geostrophic velocity, and particle path predictions. 

• MODIS SST data provided a means to monitor the position of the loop current; 
however, this method is not useful in the warmer summer months. 

• Particle path predictions on surface waters were effectively produced using the 
Runge-Kutta integration method and tri-cubic interpolation. 

• Wind data from NDBC buoys provided useful information for inferring the 
dispersal of marine debris as it nears the coast. 

• Marine debris may negatively impact critical habitat for many endangered and 
threatened coastal species, such as the Kemp’s ridley sea turtle. 

Limitations: 

• The loop current is not the only force distributing debris throughout the Gulf; 
prevailing wind patterns and long shore transport influence debris dispersal in 
nearshore environments. 

• MODIS SST was only useful during cooler months and is also hindered by cloud 
cover. 
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Future Work 



Possible future work: 

• Collect and analyze data showing marine debris accumulation along 
additional Gulf of Mexico coastlines to compare with PINS 


Use NASA data from the SeaWinds instrument onboard the 
QuickSCAT satellite for additional retrospective wind data analysis 


Include data on longshore transport rates along the Texas Gulf Coast 



Include Environmental Sensitivity Index (ESI) data for PINS to 
determine marine debris impacts to select species and habitats 

Include wind speed and direction into the Matlab codes to aid in 
predicting particle paths in a nearshore environment 
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